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We show that finite size superconductors have a spectrum of states at extremely low energy, i.e. 
inside the superconducting gap. The presence of this thin spectrum is a generic feature and related 
to the fact that in a superconductor the global phase rotation symmetry is spontaneously broken. 
For a strong coupling superconductor we find the spectrum by exactly solving the Lieb-Mattis type 
Hamiltonian onto which the problem maps. One of the physical consequences of the presence of 
thin states is that they cause quantum decoherence in superconducting qubits of finite extent. 



Introduction A most striking theorem in gauge the- 
ory was proved three decades ago by Elitzur, stating 
that local (gauge) symmetries cannot be broken spon- 
taneously The consequence of this no-go theorem 
is that the quantum mechanical expectation value of 
any local, non gauge-invariant operator must vanish. 
Naively applying the theorem to superconductors leads 
to the disconcerting conclusion that the pair amplitude 
(i/;' (r)i/)' (r)) -a manifestly non gauge- invariant quantity- 
vanishes for any position r. Taken at face value it appears 
that the picture of a superconductor as a state character- 
ized by an orderparameter with a well-defined phase and 
amplitude breaks down. Recently this observation moti- 
vated Hansson and coworkers to explore the scenario that 
superconductors are instead topologically ordered 0. 

But there is also a much more pedestrian resolution 
to the paradox that is posed by Elitzur's theorem. It 
is well known that the theorem does not forbid sponta- 
neous breaking of the global part of the symmetry. We 
exploit this caveat to construct explicitly the symmetry 
broken groundstate for a strong coupling superconduc- 
tor. This state is characterized by a well-defined global 
complex- valued orderparameter and is at the same time 
manifestly invariant under local gauge transformations. 
The appearance of the Meissner effect indicates that our 
symmetry broken state is indeed a superconductor. In 
spite of the broken phase symmetry being a continuous 
one, there are no low-energy Goldstone modes associ- 
ated with it. A well-defined energy gap develops because 
the Goldstone modes conspire with the electro-magnetic 
gauge field to create massive excitations -the well-known 
Anderson-Higgs mechanism for mass generation. 

We will prove, however, that not all states become 
gapped. A spectrum of states at extremely low energy, 
inside the superconducting gap, remains. These in-gap 
states are associated with the breaking of a global sym- 
metry in quantum systems and known in the context of 
antiferromagnets as the thin spectrum. In macroscopic 
superconductor the thin spectrum is very dense and dif- 
ficult to observe, but in finite size systems the spectrum 
displays small, a priori observable, energy gaps. Demon- 
strating their presence in any superconductor is the main 
result of this Letter. This has important physical con- 



sequences, one of which is a limit to the coherence of 
Cooper-pair box qubits. 

Below we start by constructing the symmetry broken 
state of a strong coupling superconductor, in which elec- 
trons form local pairs. The advantage of the local pairing 
Hamiltonian is that it provides a manifestly gauge invari- 
ant description of a superconductor. It also allows us to 
solve the resulting symmetry breaking Hamiltonian ex- 
actly by directly mapping it onto a modified Lieb-Mattis 
model. The fact that we consider this particular strong 
coupling model for a superconductor does not the affect 
the generality of our results on spontaneous symmetry 
breaking (SSB) and the associated thin spectrum. As 
there is no phase transition between strong and weak cou- 
pling superconductivity, these systems are identical from 
the point of view of symmetry. Or, in more specific terms: 
the thin spectrum emerges from a global symmetry break- 
ing and is therefore independent of the microscopic form 
and range of local interactions between electrons. 

Local Pairing Superconductor By means of the nega- 
tive U Hubbard Hamiltonian in presence of an electro- 
magnetic vector potential we investigate how SSB occurs 
in superconductors. This Hubbard Hamiltonian has been 
studied extensively Q and it is well known to have a su- 
perconducting groundstate. It reads 
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where ct creates an electron on site j, 8 connects neigh- 
boring sites on a cubic lattice and rij counts the num- 
ber of electrons. The amplitude of the hopping inte- 
gral is t, —\U\ the strength of the local attraction be- 
tween electrons and \i the chemical potential, which in 
general breaks particle-hole symmetry. The coupling 
of the charges to the vector potential comes about via 
the Peierls construction: while hopping from j to j + 6 
an electron picks up a phase ipj that is proportional to 
the electromagnetic vector potential integrated along the 
bond: = j. +S A s (l)dl. A gauge transformation 
amounts to simultaneously sending A — > A + VA (equiv- 
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alently: t/>j — > V| + A?+<5 ~ Aj) an d Cj — > e iA *Cj. The 
Hamiltonian is manifestly gauge invariant. We explicitly 
include the free electro-magnetic held in the Hamiltonian. 
In the next stage of the calculation this term is needed 
to generate a mass for the Goldstone modes. 

In the strong coupling limit, where U 3> t, the elec- 
trons form strongly bound pairs and for the low-energy 
dynamics of the system we can restrict ourselves to the 
lowest Hubbard sector. In that sector sites are either 
empty, or doubly occupied. Sites with single electrons 
are only virtually allowed which gives rise to pair-pair 
interactions. The effective pair dynamics can be com- 
puted from a second order perturbation expansion and is 
given in terms of pseudospin 1/2 operators as 

3,8 
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where J = 2t 2 /\U\, and = ct »ct , . The overall pair 
density is given by the z component of the total pseu- 
dospin Sf 0t , which can be varied by changing the param- 
eter h = \U\ — 2/i. Away from half filling, when h ^ 0, 
the global SU(2) symmetry of the Hamiltonian is broken 
down to a U(l) symmetry that describes the collective 
rotation of the pseudospins around the z-axis. 

Note that the prefactor of both the first and second 
term in the Hamiltonian above is J. This symmetry is 
removed when long-range Coulomb interactions between 
electron pairs are included into the model, having the 
form r V(r)SjS?j +r . For the discussion that follows, 
however, this is an inconsequential detail as such inter- 
actions preserve the global U(l) symmetry. 

The Hamiltonian above can be simplified by introduc- 
ing another set of pseudospins a that absorb the vec- 
tor potential ipj- I n this procedure, which is also used 
in weak coupling superconductivity, one defines the new 

4- —2i ib^ 4- 

pseudospins as <r| = and <rj = e ^f=a v s ST. The 
sum over j' that occurs in the definition of a + is a sum 
over the vector potentials along an arbitrary path that 
connects site j to an arbitrarily chosen origin at j = 0. 
It is easy to show that due to flux quantization this sum 
does not depend on the specific path. With this mapping 
the effective Hamiltonian for the local pairs reduces to an 
antiferromagnetic Heisenberg model in a magnetic field 

H a = J^vj ■ -h^2<j]+ H EM - (3) 
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Notice that the a pseudospin operators by themselves are 
not gauge invariant. They implicitly contain the vector 
potential and under a general local gauge transformation 
all cr-pseudospins pick up the same global phase factor 
e 2l 7fe/(°). The Hamiltonian as a whole is, of course, still 
fully gauge invariant. Following the study of spontaneous 




FIG. 1: A schematic representation of the groundstate of the 
local pairing superconductor on a square lattice. The arrows 
are semiclassical representations of the pseudospins a. 



symmetry breaking in magnets and crystals [7| we split 
up the Hamiltonian above into two parts, a finite mo- 
mentum (k) sector and a collective, zero momentum sec- 
tor. In the antiferromagnet the finite k sector contains 
the spinwave Goldstone modes which in a superconduc- 
tor become massive. These modes are disjunct from the 
zero momentum, collective sector which contains the thin 
spectrum that is needed for SSB. 

Pseudospinwaves First we consider the finite momen- 
tum (pseudo-) spinwave excitations of the Hamiltonian 
above, for which a semiclassical treatment suffices. Here 
we only give a brief account as the resulting physics is 
well known. The main point is that our Hamiltonian 
properly incorporates the Andersson-Higgs mechanism, 
generating a gap for the pseudospinwave modes @, H, @] . 

The groundstate of H a is determined by the compe- 
tition between the first and second term in the Hamil- 
tonian. The field h tends to align the spins along the 
z axis, whereas the interaction term in the Hamiltonian 
favors antiferromagnetic alignment of the pseudospins. 
The net result is, in terms of semiclassical pseudospins, 
a canted antiferromagnet in which all spins have equal 
projections on the z-axis, but are antiferromagnetically 
ordered in the xy plane (Fig. [I]). A low energy excitation 
corresponds to a long wavelength rotation of pseudospins 
around the z-axis, which amounts to a slow variation in 
the angle <f>" of each cx-pseudospin j in the xy plane. Ac- 
cording to Hamiltonian (J5]) the corresponding spinwave 
energy is proportional to J2j,s Jcos(2ip^ + (j>^ - <t>^+s), 
which in the continuum limit and for small rotations re- 
duces to an excitation energy of ^(A — V</>) 2 . The im- 
portant point is that the pseudospinwave excitations in 
the superconductor are manifestly coupled to the vector 
potential. This is made explicit by introducing a trans- 
formed vector potential A = A V</> that absorbs the 
phase rotations. Gauge invariance requires the electro- 
magnetic part of the Hamiltonian in terms of A and A to 
be identical, so that the total Hamiltonian that governs 
the elementary excitations becomes ^A 2 + Hem(A-). For 
any finite J the elementary excitations of this Hamilto- 
nian (which are combined excitations of the electromag- 
netic field and the pseudospinwaves) are massive. Via 
this Andersson-Higgs mechanism a gap is generated for 
all finite-momentum pseudospinwave excitations in our 
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superconductor. The Meissner effect is a direct physical 
consequence of this mass. 

Thin Spectrum and SSB The semiclassical approach 
above considers variations in the relative angles <f>j and 
4>j+s between neighboring pseudospins. The absolute 
angle of the pseudospins, however, is arbitrary. Conse- 
quently no unique classical groundstate exists: if all spins 
are rotated simultaneously around the z-axis by the same 
angle a different classical state results, whereas the classi- 
cal groundstate energy is invariant under such a rotation. 
In a proper quantummechanical treatment, however, the 
groundstate is both rotationally invariant and unique. In 
the following we will show that this rotational invariance 
of the quantum system can be spontaneously broken due 
to the presence of a thin spectrum. 

In exact analogy with the description of spontaneous 
symmetry breaking in for example crystals and antiferro- 
magnets 0], the collective part of the antiferromagnetic 
pseudospin Hamiltonian |3|) is given by its k = and 
k = ir part, leading to 

4j 

Hcoii = —a A ■ a B ~ ha z tot , (4) 

where we have defined two sublattices A and B between 
which the spins have antiparallel projections on the xy- 
plane. The Hamiltonian above is the Lieb-Mattis Hamil- 
tonian in presence of a uniform field h and can easily be 
diagonalized by introducing the total spin a tot = <Ja+<?b- 
The groundstate is non-degenerate and characterized by 
the total spin quantum number <j t ot and its z-projection 
a tot- As in the antiferromagnet, the groundstate does not 
break the rotational invariance of its governing Hamilto- 
nian. It is easy to see that also a set of excited states 
exists that have an energy J/N higher than the ground- 
state. These states form the thin spectrum of the sym- 
metry unbroken Hamiltonian. 

To break the symmetry an external field needs to be 
added to the Hamiltonian. Such a field is generated by 
weakly coupling our system to a second, external super- 
conductor. This leads to an effective Hamiltonian for our 
system of the form — H co u — B {a x A ~ cr%), where 

B is the effective symmetry breaking field due to the ex- 
ternal superconductor. It should be noted that in spite 
of the fact that the coupling between the our supercon- 
ductors and the external one is a gauge invariant quan- 
tity, the effective symmetry breaking field that appears 
in iJ^n is not. The great advantage that this Hamilto- 
nian can be solved exactly @, B Q counterweights the 
small complication due to the introduction of an implicit 
gauge fix. As under a gauge transformation the direction 
of the symmetry breaking field rotates around the z-axis, 
different states that are related to each other by such a 
rotation make up a gauge volume. Although they ap- 
pear to be distinct all these states in fact correspond to 
the same physical state. We will have to check that the 
thin spectrum of is made up out of distinct physical 




FIG. 2: The overlap between the thin spectrum state |n) and 
the rotated groundstate R (8) |0), as a function of the angle 
of rotation 6, for different values of n. To make this graph we 
used the values J = 10, B = h = 1 and N = 100. For higher 
values of N the graph for each n will be scaled horizontally, 
but the height of the top remains unaffected. The plot is 
symmetric under mirroring in the 6 — axis. The inset is a 
schematic representation of the dispersion relation of the low- 
energy, low-momentum states of a finite superconductor: the 
states at finite k are gapped, while the thin spectrum states 
at k — are within the gap. 



states, that are not in the gauge volume. 

The symmetry broken eigenstates of the local pair- 
ing superconductor are linear combinations of total spin 
states. The coefficients of these wavefunctions are given 
by Hermite polynomials, similar to the case of a regular 
antiferromagnet Q . The groundstate that is constructed 
in this manner has a specific absolute total phase deter- 
mined by the symmetry breaking field. Indeed it corre- 
sponds directly to the classically realized superconduct- 
ing groundstate that we considered before. The states 
corresponding to the higher order Hermite polynomials 
are extremely low in energy, and in fact collapse onto the 
groundstate in the thermodynamic limit 0, 0, fiol . [ll[ . 
They are thus easily recognized as being the thin spec- 
trum of the symmetry broken local pairing superconduc- 
tor. 

Let us now address the complication that together with 
the symmetry breaking field B we introduced an implicit 
gauge fix. One can wonder whether in this situation the 
thin spectrum states are still as physical as they are in 
the symmetry unbroken Hamiltonian ^ - the gauge fix- 
ing may have downgraded the thin spectrum states into 
the gauge volume of the superconducting groundstate. 
We can easily check that nothing of this kind happens. 
With the explicit expressions for the eigenfunctions of 
Hcofi we can evaluate the overlap between the thin spec- 
trum state |n) and the superconducting groundstate |0) 
that is rotated over an angle 9. If |n) would be in the 
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gauge volume of |0), the overlap should become unity for 
a certain rotation angle 9. As shown in figure j2| this 
overlap is equal to unity if and only if both x and 6 are 
zero. This demonstrates that an excited thin spectrum 
state is not merely a global rotation of the groundstate. 
Therefore it is not in the groundstate's gauge volume. 

BCS superconductors In weak coupling superconduc- 
tors the physical picture for symmetry breaking stays 
the same. An exactly equivalent qualitative description 
of the thin spectrum of a superconductor can be given 
within the BCS model However, it is more challeng- 
ing to obtain analytical results and closed expressions in 
this case. The physical picture comes to the fore most 
clearly if one follows Anderson by writing the standard 
BCS Hamiltonian in momentum space in terms of pseu- 
dospins [131 ] . In the superconductor the pseudospins form 
a domainwall structure around k — hp, which separates 
two ferromagnetically ordered regions with opposite pr- 
axis projections. Spontaneous symmetry breaking then 
orients the xy projection of spins in the domain wall along 
a specific direction in the plane. Again a thin spectrum 
is associated with this breaking of rotation symmetry. 

Decoherence It is known that the presence of a thin 
spectrum in mesoscopic spin qubits leads to quantum 
decoherence 0, E3]- The thin spectrum that we have 
now identified in superconductors is therefore expected 
to lead to a finite coherence time of qubits based on 
superconducting material. One such type of qubit 
that is experimentally realized is the Cooper-pair box 
qubit [3, EEL EH] • In these Cooper-pair boxes a super- 
conducting island can be brought into a superposition of 
having N and N + 1 Cooper-pairs present. Superposi- 
tions of this type can reach coherence times of up to 500 

ns E3, Ell- 
in the formalism that is outlined above it is easy to 
consider the superposition state of the superconductor 
which corresponds to the experimental one [13 ]. After 
computing the exact time evolution of such a qubit, and 
tracing over the unobservable thin spectrum states, we 
find that the coherence of the Cooper-pair box qubit de- 
cays over time. The resulting maximum coherence time 
we find to be t spon — irhN /ksT. The calculation is anal- 
ogous to the one for decoherence in antiferromagnets 0] , 
but for the superconductor N signifies the average num- 
ber of Cooper pairs on the superconducting island. Just 
as in the cases of crystals and antiferromagnets, the de- 
tails of the model (e.g. J or h) do not enter into the 
expression for the maximum coherence time, which thus 
appears as a universal timescale 0, E3] ■ Using from ex- 
periment the values N ~ 10 6 and T ~ 40 mK [II], we 
find an upper limit to the coherence time of ~ 0.5 ms for 
the experimentally realized Cooper-pair boxes. Clearly 
this timescale that is set by the presence of the thin spec- 
trum states is much larger than the timescale that is the 
current limit to coherence of the Cooper-pair boxes due 
to other environmental factors. However, it is well possi- 



ble that the limit set by the thin states will come within 
the experimental reach in the near future, either because 
the isolation from external sources of decoherence will be 
developed further, or because the size of the Cooper-pair 
box itself is reduced even more. 

Conclusions We have shown that the superconduct- 
ing state is the result of spontaneous symmetry break- 
ing. The symmetry that is broken is a global U(l) phase 
symmetry, which does not contradict Elitzur's theorem. 
In fact the resulting symmetry broken, superconduct- 
ing groundstate is still fully invariant under local gauge 
transformations. Associated with the spontaneous sym- 
metry breaking in superconductors is a thin spectrum of 
states at extremely low energies within the superconduct- 
ing gap. The presence of these low energy, global exci- 
tations leads decoherence of qubits based on supercon- 
ducting materials. The maximum time that such a qubit 
can maintain its coherence is given by t spon — TthN /ksT, 
where N is the average number of Cooper-pairs involved. 
This timescale is universal in the sense that it does not 
depend on the underlying model parameters. With the 
great technological advances in manufacturing supercon- 
ducting qubits it can be expected that this fundamental 
timescale is observed in the near future. 
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